On noise treatment in radio measurements of cosmic ray air showers 
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Precise measurements of the radio emission by cosmic ray air showers require an adequate treatment of noise. Unlike to usual 
experiments in particle physics, where noise always adds to the signal, radio noise can in principle decrease or increase the signal 
if it interferes by chance destructively or constructively. Consequently, noise cannot simply be subtracted from the signal, and 
its influence on amplitude and time measurement of radio pulses must be studied with care. First, noise has to be determined 
consistently with the definition of the radio signal which typically is the maximum field strength of the radio pulse. Second, 
the average impact of noise on radio pulse measurements at individual antennas is studied for LOPES. It is shown that a correct 
treatment of noise is especially important at low signal-to-noise ratios: noise can be the dominant source of uncertainty for pulse 
C"] height and time measurements, and it can systematically flatten the slope of lateral distributions. The presented method can also be 
>> \ transfered to other experiments in radio and acoustic detection of cosmic rays and neutrinos. 
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^ 1. Consistent definition of signal and noise 



Noise definitions applied so far in the field of cosmic ray ra- 
dio detection are originating from communication engineering. 
There, a signal usually has a power much larger than the noise, 
and lasts for a time significantly longer than its oscillation pe- 
riod. Both is not true for air shower induced radio pulses. This 
has already been investigated in the frame of self-trigger devel- 
opment ll|], where the signal-to-noise ratio plays the role of a 
threshold. For data analysis, the situation is more complex be- 
cause noise has to be defined consistently with the definition 
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of the radio pulse height, which is the maximum of the field 
strength, in the case of LOPES J2D. 

Independently of the specific signal and noise definitions, the 
following consistency criterion is demanded: 



for true signal — 



measured signal 



= 1 



noise 



(D 



The consistency criterion is supposed to hold only on aver- 
age, because the noise level at the signal time can by chance be 
larger or smaller than the average noise level. In addition, even 
for a positive true signal, the measured signal-to-noise ratio can 
in some cases be smaller than one, since noise can interfere con- 
structively or destructively with the air shower radio emission, 
and increase or decrease the measured signal compared to the 
true signal. 
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Figure 1: Typical noise measured with a LOPES antenna (left), and a test pulse from a pulse generator (right): sampled data points, the up-sampled trace and a 
Hilbert envelope of the trace are shown in both cases. The noise level is calculated as the weighted average of the local maxima of the envelope. This corresponds 
to the average level of the plotted step function with a step exactly in the middle between two local maxima. 
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Figure 2: Signal-to-noise ratios of pure noise for a selection of LOPES events 
without signal. 



For LOPES, a consistent definition of signal and noise has 
been found for measurements at individual antennas, e.g., to re- 
construct the lateral distribution [3]. The signal is defined as 
the amplitude (field strength) of the radio pulse which is deter- 
mined as the local maximum of an Hilbert envelope closest to 
the pulse time known from a preceding interferometric cross- 
correlation beam analysis (c.f. A). The noise level is defined 
as average amplitude in a time window (10 |xs) before the radio 
pulse, and is calculated by the mean of all local maxima of the 
envelope. Because lower local maxima are more likely to have 
a smaller distance to neighboring maxima than higher maxima, 
it is necessary to weight each maximum with the distance to its 
neighbors when averaging (fig.[TJ left). 

It has been tested that these definitions of signal and noise do 
indeed fulfill the consistency criterion. With a selection of 200 
LOPES events without radio pulse, a large sample of 120,000 
noise intervals of 10 ^is width, each, has been obtained: the in- 
tervals are non-overlapping, cover different days and times of 
the day, as well as different antennas. The average signal-to- 



noise ratio of these intervals is compatible with 1, as required 
(fig.0. 

With other definitions of noise, like the RMS of the field 
strength or its square (power), the mean of the absolute field 
strength or an unweighted mean of the local maxima of the en- 
velope, the consistency criterion is not fulfilled. However, the 
ratio between the noise levels determined by different methods 
is constant within a few percent. Thus, results obtained with a 
different noise definition could be scaled to a consistent defini- 
tion when accepting a small systematic error. 

2. Influence of noise on pulse height measurements 

The impact of noise on measurements of the pulse ampli- 
tude at individual antennas has been studied for LOPES with 
test pulses (fig. [T] right) of different width, and noise from real 
measurements. Therefore, the test pulses have been scaled with 
the LOPES analysis software to a certain amplitude A^ rue , and 
added to the noise intervals presented in the previous section. 
Afterwards, the measured signal height Ameas can be obtained 
for each pulse, yielding a relation between the average true am- 
plitude At me and the measured amplitude Ameas- To simplify 
the relation, all amplitudes have been normalized to the noise 
level, i.e., the noise level corresponds to A = 1, and Ameas is 
the measured signal-to-noise ratio. 

Because the real probability distribution of A( rue of air 
shower induced radio pulses is unknown, scaling factors for the 
test pulse heights have been generated for a flat distribution of 
A]- me . As cross-check, also an exponentially decaying distribu- 
tion has been tried, but the effect on the results is negligible. 

To correct measured pulse amplitudes for the noise influence, 
the function At me (Ameas) is required, which has been obtained 
by the following procedure. The test pulse data, which con- 
sist of 120,000 samples with known Atr U e and corresponding 
Ameas, have been sorted into bins. Thereby, each bin cov- 
ers a certain interval of Ameas- The mean A true of each bin 



2 



o 

"D 




LOPES data 
amplitude parametrization 



4 6 8 10 12 
measured signal-to-noise ratio 



test pulse data * 

parametrization 

impact of noise correction « 
calibration uncertainty 




2 4 6 8 10 12 14 

measured signal-to-noise ratio 



Figure 3: Relation between the true signal height Atme and the measured 
signal-to-noise ratio A me as- The error bars correspond to the standard devi- 
ations of the binned data, i.e. AAt me , and not the the uncertainty of the plotted 
mean. 



is then the average true amplitude corresponding to the mea- 
sured amplitude Ameas of the bin. At the same time, the stan- 
dard deviation AA^g of each bin can be taken as error estima- 
tion of the true amplitude (fig. [3). Other methods to determine 
A t rue<Ameas) failed. The inverse function of Ameas(A t rue)> 
which would be available directly, is not defined for Ameas < 1 ■ 
Using confidence intervals instead of mean and standard devi- 
ation yields meaningless results for Ameas close to 0. More 
details about the method will be available in |5[ . 

The validity of the relation At me (Ameas) has been checked 
for various systematic effects, and no significant dependencies 
could be found. In detail, the following effects have been stud- 
ied: the (up-)sampling rate, the shape of the test pulse, the an- 
tenna type and polarization used for noise measurement. Sum- 
marizing, any possible effects are negligible against the size of 
the error AAt rue , and against the calibration uncertainty due to 
environmental effects (~ 5 %, c.f. |0]). 

The following parametrization has been found for Aj- me : 



^true - yj- 



A meas ~ 1 for A meas > 2 



and for low signal-to-noise ratios: 
^true - a + b ■ A me£ 



for Ameas < 2 



(2) 



(3) 



with a = 0.4628 + 0.0066, b = 0.2491 +0.0092, and c = 2.349 + 
0.048, determined by a fit, with a forced connection to the first 
formula at Ameas = 2. 

This means that at high signal-to-noise ratios, the power of 
the measured signal is in average the sum of the power of the 
noise and the power of the radio pulse, i.e., noise generally is 
more likely to increase the signal amplitude than to decrease 
it. However, at low signal-to-noise ratios the behavior is not 
trivial. This demonstrates that a detailed study of the noise in- 
fluence, like performed here for LOPES, is indeed necessary, 
especially because most LOPES events contain at least some 
antennas with signal-to-noise ratios < 2. 



Figure 4: Relative uncertainty of the true pulse height ( AAt me /At rue , x), in 
each bin of the measured signal-to-noise ratio A me as ■ This is compared with 
the error one would make, if the measured pulse amplitude would not be cor- 
rected for the noise influence ((A[ rue -Ameas)Mtrue> +)• and w i m me calibra- 
tion uncertainty due to environmental effects (horizontal lines). 
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Figure 5: Example lateral distribution from a typical selection with and without 
correction of the noise influence. 



The noise dependent, statistical error of Ameas which is de- 
termined as standard deviation AAt rue , like explained above, 
clearly exceeds the calibration uncertainty for signal-to-noise 
ratios < 10 (fig.|U). Like for the amplitude itself, also the error 
AAmeas is parametrized differently for low and high signal-to- 
noise ratios: 

A^true = d + e ■ exp (+A m eas) for A m eas < 1 -68 (4) 
A^true =f + g- exp (-Ameas) for A m eas > 1-68 (5) 

with d = 0.3103 ± 0.0096, e = 0.0647 + 0.0029, / = 0.6162 + 
0.0010, and g = 0.213 ± 0.018, obtained from a fit. 

3. Influence on lateral distributions 

Since the amount by which noise in average increases or de- 
creases the amplitude of radio pulses depends on the signal-to- 
noise ratio, it also has an impact on the slope of lateral distri- 
butions of the air shower induced radio emission. For a typical 
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Figure 6: Average deviation between the measured time of the envelope maxi- 
mum after adding noise, and the original time before adding noise, and the time 
calibration uncertainty as reference. 

selection of LOPES events (e.g., the one used in [3]), noise sig- 
nificantly flattens lateral distributions (fig.[3]l, except for events 
with very high signal-to-noise ratios at all antennas. 

If an exponential function A{R) = eo • exp(/?//?o) is fitted to 
each lateral distribution with a slope parameter Rq, noise typi- 
cally increases Rq by 10 - 20 %, which is not negligible com- 
pared to other uncertainties. The effect on the amplitude param- 
eter is smaller and vanishes only for events with high signal-to- 
noise ratio at all antennas. Fortunately, it is now possible to 
correct for the influence of noise in every single measurement, 
at each individual antenna, by the parametrization formulas pre- 
sented above, to obtain the 'true' lateral distributions. 

4. Influence on pulse arrival time measurements 

The influence of noise on the measurement of pulse arrival 
times has been studied in a similar way like the influence on 
pulse amplitudes, defining the pulse arrival time at an antenna 
as the time when the pulse amplitude is measured. Noise ran- 
domly shifts the true pulse arrival time of test pulses to earlier 
or later times, and no significant tendency to either one is ob- 
served. 

The mean of the absolute shifts of the pulse time at a cer- 
tain signal-to-noise ratio is defined as the noise dependent time 
uncertainty At (fig.|6]l. For signal-to-noise ratios < 10, it can 
be much larger than the time calibration uncertainty, which is 
about 2 ns for pulse arrival time measurements at individual an- 
tennas |0. For the interferometric cross-correlation beam ana- 
lysis used in LOPES, the impact of noise will be smaller, be- 
cause shifts at different antennas will average out, but a quanti- 
tative study is beyond the scope of this paper. 

Several possible systematic effects have been examined. As 
expected, At does neither depend on the antenna type nor its po- 
larization. It does depend on the sampling frequency, but the ef- 
fect becomes negligible against the time calibration uncertainty 
of LOPES, if data are up-sampled to at least 640 MHz. Unfor- 
tunately, At depends on the pulse shape, but no correlation with 



easy accessible parameters like pulse height or width could be 
found. For this reason, and since it is unknown which test pulse 
shape does best describe the real cosmic ray induced pulses, the 
following parametrization of At takes into account the average 
behavior of all tested pulse shapes: 

A? = A m ^ s • 20.5 ns for A me as > 1 -8 (6) 

whereby Ameas is the measured signal-to-noise ratio, and all 
parameters have been adjusted by hand to fit the data. The ar- 
rival time uncertainty At does not become arbitrary large for low 
signal amplitudes, because the time interval for signal search 
depends on the preceding cross-correlation beam-forming. For 
high signal-to-noise ratios, the behavior is consistent with expe- 
rience from beam-forming where the resolution improves with 
the signal-to-noise ratio. 

5. Conclusions 

Treating noise correctly in measurements of radio pulses 
emitted by air showers is especially tricky at low signal am- 
plitudes. Nevertheless, it is important because events close to 
the detection threshold will always contain antennas with low 
signal-to-noise ratios. For instance, it has been shown that 
noise systematically flattens lateral distributions measured with 
LOPES. 

Alternatively, only events at high signal-to-noise ratios could 
be studied, where the noise influence becomes negligible 
against calibration uncertainties. Then, the signal could be de- 
fined as the integrated power of the radio pulse, if the integra- 
tion time is large against the pulse width. This is a tempting 
approach, because this method could be realized directly in ana- 
log electronics, and would allow cheaper and easier designs of 
radio detectors. However, this will be paid with a higher de- 
tection threshold, and is no option for current experiments like 
LOPES or AERA|H. 

Independently of the experimental design, attention must be 
paid to define signal and noise consistently in any analyses. 
Consistent definitions have been presented for LOPES for pulse 
measurements at individual antennas. A corresponding study of 
the noise influence on beam-forming analyses is more complex, 
and might better be performed with simulated pulses. Neverthe- 
less, the presented method and results for the noise impact at in- 
dividual antennas can probably be transfered to any experiment 
based on radio or acoustic arrays, where the signal consists of a 
short, bandwidth limited pulse, and noise can interfere in both 
ways, destructively and constructively with the signal. 
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